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Abstract 



Recent Atomic Force Microscope (AFM) nanoindentation experiments measur- 
ing mechanical response of the protein shells of viruses have provided a quantita- 
tive description of their strength and elasticity. To better understand and interpret 
these measurements, and to elucidate the underlying mechanisms, this paper adopts 
a course-grained modeling approach within the framework of three-dimensional 
nonlinear continuum elasticity. Homogeneous, isotropic, elastic, thick shell mod- 
els are proposed for two capsids: the spherical Cowpea Chlorotic Mottle Virus 
(CCMV), and the ellipsocylindrical bacteriophage 4>29. As analyzed by the finite 
element method, these models enable parametric characterization of the effects of 
AFM tip geometry, capsid dimensions, and capsid constitutive descriptions. The 
generally nonlinear force response of capsids to indentation is shown to be in- 
sensitive to constitutive details, and greatly influenced by geometry. Nonlinear 
stiffening and softening of the force response is dependent on the AFM tip dimen- 
sions and shell thickness. Fits of the models capture the roughly linear behavior 
observed in experimental measurements and result in estimates of Young's moduli 
of ^280-360 MPa for CCMV and «4.5 GPa for 029. 

Key words: viruses; atomic force microscopy; computer simulation; continuum 
mechanics; elasticity 



Introduction 



As confirmed by a scan of the RCSB Protein Data Bank ( i), the structural biology 
techniques of X-ray crystallography and cryo-electron microscopy (cryo-EM) have 
enabled determination of the shape and structure of a multitude of macromolecules, 
in many cases locating individual atoms with resolution up to a few angstroms. In 
particular, these methods have been effective in elucidating the structures of the 
protein shells (capsids) of many viruses (|2j). However, despite (or perhaps because 
of) the wealth of detailed structural information made available by these methods, 
it remains a significant challenge to understand and predict the overall structural 
mechanics of capsids based on models of fundamental atomic interactions. Global 
deformations of capsids, such as those implicated in viral maturation (3) and infec- 
tion processes (0,0), can entail large-scale coordinated motions in which all atoms 
are involved. The shear number of atomic degrees of freedom in a viral shell makes 
molecular dynamics simulation of such processes unreasonable. In order to under- 
stand coarser-scale capsid mechanics, experimental and theoretical approaches are 
needed that describe the capsid structure in reduced terms. 

Atomic force microscopy (AFM) is such a tool, able to explicate virus me- 
chanics by measuring the force response due to indentation of a capsid. Recent 
AFM studies probing the mechanics of viral capsids da |7|) have shown that they 
can be strong and yet highly elastic, even under significantly large deformations. 
A particularly interesting hallmark of these experiments is the linearity of force- 
indentation response. Prompted by this observation, Ivanovska et al. interpreted 
AFM experiments on the cp29 bacteriophage by building theoretical models based 
on (linearized) small-strain, thin-shell, continuum elasticity, thus enabling quanti- 
tative estimation of three-dimensional elastic constants (e.g., the Young's modulus) 
for the capsid. More recent experiments have measured a similarly linear force- 
indentation response for a completely different virus, Cowpea Chlorotic Mottle 
Virus (CCMV) Q), apparently affirming the appropriateness of linearized elasticity 
modeling. However, though linear-response theory does offer a description of these 
AFM experiments, it is not at all clear why a linear large-deformation response is 
observed in the first place. Indeed, in traditional solid mechanics, large deforma- 
tions are frequently accompanied by a nonlinear force response, even when the 
material is perfectly elastic. The nonlinear elasticity of rubber materials is perhaps 
the most obvious example <H). 

In this paper the a priori assumption of linearity is abandoned, and more gen- 
eral elastic continuum capsid models are developed and studied. Capsids are mod- 
eled in the context of finite-deformation hyperelasticity, wherein strains are not as- 
sumed to be small, such that the nonlinear effects of large displacements, rotations, 
and strains are considered. Previous linearized elasticity models of capsid indenta- 
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tion (0, 0) have also relied on reduced thin-shell theory (fiol 11 ), under which the 
three-dimensional equations of elasticity are reduced to two dimensions, facilitat- 
ing analytical solutions for simple loadings and geometries. However, thin-shell 
theory, as its name implies, is applicable only for shell structures with a thick- 
ness that is much smaller (by at least a few orders of magnitude) than the overall 
structural dimensions. In particular, thin-shell theory precludes transverse shear 
deformation, which though negligible for thin shells, is an important factor for 
thick shells. The nominal thickness of CCMV is just over 10% of its outer diame- 
ter, likely putting it outside the range of applicability for thin-shell analysis. Two 
options remain for analysis of a thick shell like the CCMV capsid. The first is to 
resort to shell theories which allow for shear deformation (10), and the second is 
to work within the general framework of 3-D continuum elasticity. In the present 
work, focus is on the latter, more general approach, and its differences with the 
former are briefly examined . Both of these strategies are difficult to employ in a 
closed-form analytical manner, and typically for scenarios involving large defor- 
mations and complex geometries or boundary conditions they require numerical 
solutions. The capsid indentation experiment is one such scenario, involving both 
large strains and geometric nonlinearities due to contact between the capsid and the 
substrate/AFM. In this work numerical solution is performed by the finite element 
method dl2l) . the most popular and robust analysis technique for elasticity problems 
in engineering. Calculations are carried out using ABAQUS Jl3l) . a commercially 
available nonlinear finite element analysis software package. 

Finite element simulations of AFM nanoindentation experiments provide a 
means for systematic study of the effects of model parameters on computed force- 
indentation response curves, thereby enabling validation and comparison to ex- 
periments. The inputs to the continuum elasticity models are of constitutive and 
geometric natures. To understand the influence of constitutive modeling choices, 
the classical (Hookean) linear stress-strain law is compared with nonlinear hypere- 
lastic constitutive laws. Also, the dependence of force-indentation curves on capsid 
dimensions and AFM-tip loading geometry is demonstrated. The results show that 
the nonlinear contact interaction of the AFM tip with the capsid has important 
qualitative and quantitative influence on the structural response. 

The experimental studies (la, 13) that have motivated this work measured the 
nanoindentation response of two viruses: Cowpea Chlorotic Mottle Virus (CCMV) 
and bacteriophage 029. CCMV is a roughly spherical plant virus of about 28 nm 
in diameter. Its protein capsid, which protects a single-stranded RNA genome, is 
assembled from 180 identical protein subunits arranged into a truncated icosahe- 
dron with a triangulation number T=3, according to the Caspar-Klug classifica- 
tion scheme for icosahedral viruses (14). 029 is a bacteriophage (i.e., a virus that 
attacks bacteria) with a roughly ellipsocylindrical capsid that protects a double- 
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stranded DNA genome. The 029 capsid is comprised of 235 protein subunits, 
arranged into a prolate shell with a center cylindrical region and two icosahedral 
end caps 

The following section describes the mechanics framework and the finite ele- 
ment formulation for the models, and presents the idealized structural models of 
CCMV and 029. Next, results are presented of simulations of nanoinentation of 
CCMV along with some estimates of effective elastic moduli. Because its more 
symmetric shape makes it a simpler candidate for modeling, CCMV is focused on 
as the subject of a series of parametric studies on model parameters. First the sensi- 
tivity of the model to changes in constitutive laws is assessed. Despite the presence 
of large strains in simulations, no significant difference in overall force-indentation 
response is observed for the three different models. Secondly the influence of the 
AFM tip geometry and shell thickness on the overall response is examined. These 
geometric properties are shown to produce both stiffening and softening nonlin- 
earities in the structural force response. For comparison purposes, the computed 
indentation response of a 029 is shown, and an estimate of its effective elastic 
moduli is made. 

Models and Methods 
Modeling capsid geometry 

The native atomic structure of CCMV was determined to within 3.2 A by Speir 
et al. (4) using X-ray crystallography and cryo-electron microscopy. In this study, 
the capsid was shown to be a T = 3 icosahedral shell, composed of protein subunits 
organized with regions of five-fold symmetry (pentamers) and six-fold symmetry 
(hexamers) that form protrusions from the capsid surface, making the thickness 
of the shell highly nonuniform. For the outer radius of CCMV, Speir et al. (0) 
measured a maximum value of RS^L = 14.2 nm and a minimum value of = 

max ■ ■ min 

12.0 nm. The average outer Radius is R™ g = 13.2 nm did) . For the inner surface 
of the capsid they specify a minimum radius of R^ in = 9.5 nm and an average 
radius of R™ g = 10.4 nm. Here CCMV is idealized as a uniformly thick spherical 
shell as shown in Fig.^a). It is not immediately clear what uniform dimensions for 
the idealized shell would best represent the highly nonuniform physical dimensions 
of the actual capsid. As a starting point, average radius values are adopted for the 
inner and outer surfaces, which from (|4|) are R m = 10.4 nm and i? out = 13.2 nm, 
such that the idealized thickness is t = 2.8 nm. Given the arbitrariness of adopting 
average dimensions for the idealized model, a parametric study is done to examine 
the influence of these choices on the model results. 
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The second virus of the study, bacteriophage 029, has a less symmetric but 
more uniform structure, as determined by Tao et al. (1151) using cryo-electron mi- 
croscopy. Its capsid is formed by two T=3 icosahedral caps connected by a cylin- 
drical region of hexamers. The average thickness of the capsid is 1.6 nm. The two 
outer diameters are known to be 54 nm from the apex of one endcap to another and 
42 nm across the cylindrical region The 029 capsid is modeled by a center 
cylindrical region 24 nm in length connecting two ellipsoidal endcaps. Ellipsoidal 
endcaps (rather than spherical) are chosen as there is a disparity between the mea- 
sured dimension from the center to the cylindrical section (21 nm) and from the 
center to the apex of the endcap (15 nm). If the endcaps were made spherical, and 
the length of the cylindrical region chosen to match the reported 24 nm value, then 
the apex to apex diameter of the model would be much larger than the reported 
value. Model dimensions are summarized in the schematic diagrams in Fig.Q 




(a) CCMV (b) 029 



Figure 1: Schematic of idealized capsid models. Dimensions are averaged from 
the results of structural studies of CCMV (0) and 029 (Q). 

Theoretical background and computational methodology 

Following the standard definitions of finite deformation continuum mechanics 
HiL the deformation of a body is described by the one-to-one deformation map- 
ping ip : Bq — » M 3 , which maps material point at position X G Bo in the reference 
configuration of the body to point x = ip(X) in the deformed configuration of 
the body. The deformation gradient tensor is denoted F = Vip, the right Cauchy- 
Green deformation tensor C = F T F and the Green strain tensor E = |(C — 1). 
The volume change ratio of current to reference volume is given by the determi- 
nant of the the deformation gradient, J = 4ft- = det F. Also, in the context of 
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compressible large-strain elasticity, it is helpful to introduce isochoric (volume- 
preserving) deformation tensors F = J~ 1//3 F and C = F T F. 

Constitutive theory for hyperelastic materials postulates the existence of a strain 
energy density function W, which gives the stored elastic energy per unit refer- 
ence volume at every material point in the body. The stress response of a hy- 
perelastic material is then given in terms of the first Piola-Kirchhoff stress tensor 
P = dWj dF or the second Piola-Kirchhoff stress tensor S = dW/ <9E. 

Mechanical equilibrium can be enforced by minimizaton of the total mechani- 
cal energy (including work of contact forces). The finite element method approx- 
imated the nonlinear elasticity problem by interpolating the deformation mapping 
among the deformed positions of the vertices (nodes) of a polyhedral mesh, and 
solving for the nodal positions which best minimize energy. The finite element 
equilibrium equations are nonlinear in the unknown deformed nodal positions x a 
because of the nonlinearities in the constitutive expressions for stress and (for the 
present context) in the contact conditions which define external forces. 

Constitutive laws 

In this work, three particular hyperelastic models are employed: St. Venant-Kirchhoff, 
Neo-Hookean, and Mooney-Rivlin. The so-called St. Venant-Kirchhoff model, is 
given by the strain energy function 

^ = ^(trE) 2 + /[ / (trE 2 ). (1) 

The chief appeal of this material model is that the resulting (second Piola) stress 
response is linear in the Green strain, S = Ao (tr E) 1+2/^oE, and is, in a sense, akin 
to Hooke's law for isotropic (infinitesimal-strain) linear elasticity. Linearization 
of this stress response for vanishingly small values of the displacement gradient 
H = F — 1 identifies Ao and as the initial Lame constants. 

For isotropic compressible hyperelastic materials it can be convenient to define 
the strain energy function in a decoupled representation as 

W = W vol {J) + W iso {I 1 ,I 2 ), 

where W vo \(J) describes the volumetric response, and Wiso(h-h) describes the 
isochoric response in terms of the first two principal invariants of C. The simplest 
example of such a constitutive law, the Neo-Hookean model defines the isochoric 
response to be simply linear in the first invariant: 

W h0 = ff(h-3). (2) 
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Emerging from the statistical mechanics of a Gaussian chain , the Neo-Hookean 
form has some justification as a highly idealized constitutive model for polymers 
and proteins. Linearization of the Neo-Hookean stress response yields hq as the 
initial shear modulus. The Mooney-Rivlin model can be though of as an exten- 
sion of the Neo-Hookean to include the lowest-order dependence on the second 
invariant: 

W iso = Cio {h - 3) + Coi (J a - 3) . (3) 
Linearization of Mooney-Rivlin for small strains gives an initial shear modulus 
of no = 2(Cio + Coi)- F° r both Neo-Hookean and Mooney-Rivlin models, we 
choose a simple quadratic volumetric response W m \ = k${J — l) 2 /2, such that 
linearization renders ko = Ao + §/xo as the initial bulk modulus, consistent with 
initial Lame constants, Ao and fiQ. 

Simulation of capsid nanoindentation 

AFM nanoindentation of capsids is simulated by the finite element method, using 
the commercially available package ABAQUS dl3h . The models described above 
are meshed with finite elements, and the capsid meshes are compressed between 
substrate and AFM tip. The substrate and AFM tip are both modeled as rigid. The 
substrate has the geometry of a flat plate, whereas the AFM tip is given a spherical 
shape with a nominal radius of about 14 nm, consistent with geometry of the tips 
used in the experiments of Ivanovska et al. \6) and Michel et al. Q). Though the 
spherical CCMV model obviously has no intrinsic orientation, the same is not true 
for the ellipsocylindrical </>29. Again motivated by the experiments of Ivanovska 
et al. (@), the (f>29 model is oriented such that compression is simulated along the 
short axis of the capsid. Noting presence of symmetries in the capsid and loading 
geometries, it seems reasonable to expect deformation to obey the same symme- 
tries. Hence only one-quarter of each capsid is meshed, and symmetry boundary 
conditions along the boundaries are enforced. Meshed model assemblies are shown 
in Fig. El 

Contact of the rigid tip and substrate with the deformable capsid is modeled as 
frictionless, such that Lagrange-mulitiplier contact forces are introduced normal to 
the contact interfaces, in order to prevent interpenetration. The progressive com- 
pression of the capsid is achieved by an incremental, quasi-static, displacement- 
controlled process, where each increment involves a small upward displacement 
of the rigid substrate followed by iterative solution of the nonlinear finite element 
equilibrium equations consistent with contact constraints. For each incremental 
displacement of the substrate, the total contact force at each contact interface is 
computed, allowing construction of a force-deflection curve for the indentation 
process. 
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Figure 2: Finite element meshes of CCMV (left) and <p29 (right) Capsids. Dimen- 
sions of the capsids are shown in Fig. [2 AFM tip is modeled as a rigid hemispher- 
ical shell with R = 14 nm, substrate as a rigid flat plate. (Capsids are shown not 
on the same scale.) 

Nanoindentation of CCMV 

As described in the previous section, parameters defining the spherical model of 
CCMV are of two types: geometric and constitutive. The simulations described 
here are designed to separately and quantitatively understand the effects of these 
parameters on capsid indentation. The influence of capsid dimensions is considered 
by fixing the radius of the midsurface of the shell (equidistant from the inner and 
outer surfaces) at R = 11.8nm, while varying the thickness from 1 nm to 5 nm 
(roughly reflecting the range of physical thickness over the actual capsid shell). A 
range of AFM tip geometries is also considered, from a rather sharp tip of radius 
3.5 nm to a completely flat tip (i.e., infinite radius). 

The three consititutive theories described above — St. Venant-Kirchhoff, neo- 
Hookean, and Mooney-Rivlin — are parameterized in a manner such that they lin- 
earize consistently with initial Lame constants related to initial Poisson ratio v and 
initial Young's modulus E by 

_ vE E 

Informal comparisons showed that the model results were quite insensitive to the 
choice of Poisson ratio. Because proteins tend toward behaving incompressibly in 
their elasticity, an only slightly compressible Poisson ratio of v = 0.4 is chosen 
here for all simulations. As can be seen from the definitions earlier, because the 
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(a) J "in 



(b) 7nm 




(e) 9 nm 



Id) 12 nm 



Figure 3: Indentation of spherical model of CCMV with dimensions from Fig. Q 
and a constitutive response defined by the St. Venant-Kirchhoff theory with E = 
250 MPa, and v = 0.4. Deformed capsid shapes are shown for several values of 
substrate displacement. Color contours indicate the von Mises stress. 



simulation is displacement-controlled, the Young's modulus acts effectively as a 
simple proportionality factor for the stress response of the material (and hence 
also for the force response of the entire shell). Therefore, simulation results can 
be normalized by E, avoiding the need for running a series of simulations over a 
range of values. 

Parametric studies of geometric and constitutive effects are described in the 
following section. As a representative example, consider a simulation of CCMV 
with the nominal dimensions from Fig.^(outer radius of 13.2 nm and a thickness of 
2.8 nm), described by St. Venant-Kirchhoff elasticity, indented by a spherical rigid 
AFM tip of radius 14 nm. For a capsid with Young's modulus of E = 250MPa 
(on the order of what's expected for typical protein material (20)), Fig. |3 shows 
a sequence of deformed capsid shapes, color coded by von Mises stress (EH) at 
several intervals during the indentation process. 

Fig-IUshows that as AFM indentation increases during simulation, elastic stresses 
build and are distributed through the thickness of the capsid wall (as indicated by 
the contours of von Mises stress). At the outset, the stress is highest at and around 
the point of contact between the AFM tip and the top of the capsid, and is more 
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uniformly distributed through the thickness away from the contact region. As the 
substrate displacement continues, the stress is increasingly concentrated on the in- 
side surface of the capsid. In general at later stages in the indentation, the von 
Mises stress decreases in magnitude through the thickness, taking its lowest value 
at the center of the capsid wall, consistent with the interpretation that bending be- 
comes a more dominant load-carrying mechanism. For this and all other CCMV 
simulations, indentation is carried forward until the substrate displacement reaches 
30% of the capsid diameter, approximately the point of failure in experiments O). 
At this point in this representative simulation, the largest von Mises stress is about 
90 MPa and the maximum magnitude of the principal strains is about 15%, clearly 
outside the range of applicability of small strain linear elasticity. Since many mate- 
rials have ultimate strengths that are in the range of 1-10% of their Young's moduli 
the maximum stress of 90 Mpa far exceeds this rule of thumb for indicating 
failure. This suggests the capsid is qualitatively similar to rubber materials, which 
typically have higher ultimate strengths and lower Young's moduli. 

Given the magnitude of strains, along with the significant changes in contact 
geometry throughout the simulation depicted in Fig. |3] it is somewhat surprising 
that the corresponding force-indentation curve, shown in Fig. 01 exhibits only very 
subtle nonlinearity. For very small displacements of the substrate, there is a notice- 
able positive curvature in the force response. This corresponds to the initial part 
of the indentation, where deformation gradually spreads through the thickness of 
the shell as the contact area between tip and capsid grows. In the limit of small 
deformations, this should lead to a Hertz-like contact force scaling, which from 
linear elasticity is expected to take the form F ~ d 3 ^ 2 did I 111) , where F is the to- 
tal contact force and d is the indentation or change in capsid height. This expected 
scaling is consistent with initial convexity of the simulated force response. 

At larger indentations, the curve becomes remarkably linear over a large range, 
before eventually softening around d 7 — 8 nm as bending deformations of the 
shell become more severe. It is in this range in the simulation that the contact 
pressure near the center of the AFM tip contact region begins to decrease to the 
point where eventually full contact is lost and a gap opens between the AFM tip 
and the apex of the capsid. The contact region is from then on defined by a circular 
ring. Fig. |3 a) shows a close-up side view of the AFM tip and the top of the capsid 
illustrating the separation of the capsid away from the AFM tip. Fig. |3b) shows 
the contact pressure between the AFM tip and capsid surfaces at indentations of 
d 8 nm and d ~ 9 nm, just after the capsid and AFM tip separate. The blue 
area in the center of the top portion of the capsid represents zero contact pressure, 
and the circular ring area where the AFM tip and capsid are still in contact can 
clearly be seen. Zero contact pressure is indicative of a loss of contact several 
nanometers before the separation (seen in Fig.|5t) is visible. The inner gray region 
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Figure 4: Force response for indentation of spherical CCMV model from Fig.|3] 



in Fig. |5la) is the area of the two surfaces that are no longer in contact. This sepa- 
ration of the capsid apex from the AFM tip geometrically resembles the so-called 
"snap-through" buckling of arches and thin shells did l22h . This nonlinear contact 



effect is not to be confused with a geometric instability or equilibrium bifurcation, 
which would be marked by a drop in the force-indentation curve. However, the 
softening of the force-indentation curve is qualitatively similar to "snap-through" 
behavior, and consistent with a transition of the primary load-bearing mechanism 
from stretching to bending of the shell surface. The absence of a drop in force for 
this reversibly elastic model perhaps suggests that the failure experimentally ob- 
served by Michel et al. (Q) may be related to instability in the constitutive response 
of the protein material. 



Parametric studies 

Constitutive laws. As shown in the previous section, AFM indentation is capable 
of inducing significant strains in a capsid shell. The experimentally-demonstrated 
fact that CCMV capsids remain elastic despite such deformations, puts them in 
a very special class of materials. Phenomenological constitutive theory for large- 
strain nonlinear-elastic materials has traditionally focused on the response of rubber- 
like materials. Two of the lower-order theories for rubber-like materials: neo- 
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(a) Visible loss of contact. d~ 12 nm 



(b) Contact pressure 



Figure 5: At higher deformations, the apex of the CCMV capsid is observed to 
separate from the indenting AFM tip. Separation is preceded by a drop in contact 
pressure near the apex. 
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Hookean and Mooney-Rivlin, are employed here to get a sense for how impor- 
tant constitutive details are relative to geometric nonlinearities in determining the 
overall shape of the force-indentation curve. Fig.[6]compares force-indentation re- 
sponses for capsids modeled with neo-Hookean and Mooney-Rivlin to that of the 
linear St. Venant-Kirchhoff model. The parameters for all of these simulations are 
identical to those the representative CCMV model above, with the dimensions of 
the shell taken from Fig.Q a Poisson's ratio of v = 0.4 and an initial Young's mod- 
ulus of E = 250MPa. The coefficients of three-parameter Mooney-Rivlin model 
are determined by specifying the fraction in addition to E and v. Three values 
of this fraction are considered, such that the first invariant is weighted respectively 
less, evenly, and more than the second invariant. Fig. |6] clearly demonstrates that 
the force-indentation curves for all of these constitutive models are practically co- 
incident throughout the entire deformation. This clarifies that the details of the con- 
stitutive response are not centrally important to the shape of the force-indentation 
curve. It suggests rather that the geometry of the shell, and the geometry of loading 
are more influential in determining the indentation response. Furthermore, these re- 
sults support the appropriateness of the linear stress-strain law. Apparently, despite 
the presence of rather large strains, constitutive nonlinearities are masked from the 
overall structural response. This is likely due to the fact that the strains over much 
of the capsid are significantly smaller than the maximum value, such that locally 
even the nonlinear models are in the linear response regime. 

AFM tip size. In the experiments of (0,0) capsids rest on a flat substrate and are 
compressed by an AFM tip which is roughly hemispherical. The AFM tip radius 
is on the order of the outer radius of the CCMV capsid, but exact determination of 
the dimensions is difficult. Hence there is a need to determine the importance of 
the modeled AFM tip size on the resulting contact behavior. Simulations were per- 
formed varying the AFM tip radius from one-quarter of its nominal size (14 nm) 
to infinite (treating the AFM tip as a second flat plate), and the resulting contact 
force curves are compared in Fig. Q At displacements of 7 — 8 nm, when failure 
is experimentally observed, the contact force varies by rs 15%. When approach- 
ing either extreme in the size of the AFM tip, be it of infinite radius (a flat plate) 
or nearing a point load (one quarter of the initial radius), the results change only 
slightly. Also notable is the subtle change in shape of the force-indentation curve 
as the tip radius is changed. Softening of the force response, noted earlier to occur 
along with a separation of the capsid apex from the tip, is accentuated by smaller 
tips which induce more bending deformation in the capsid. Indeed, over the range 
of deformation simulated, this softening is completely absent when the tip is mod- 
eled as a flat plate. In this case bending deformations near the apices is limited and 
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AFM Plate Displacement (nm) 



Figure 6: Sensitivity to constitutive law. Contact force curves for St. Venant- 
Kirchhoff, compressible Neo-Hookean, and compressible Mooney-Rivlin models. 
Representative capsid dimensions from Fig. [2 are used. 
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Figure 7: Contact force comparison for several AFM tip sizes. Representative 
capsid dimensions from Fig.[T]are used. 



the slope of the force-indentation curve is monotonically increasing. 

Capsid thickness. As noted earlier, the Young's modulus, E, acts as a simple 
proportionality factor for stress and force responses to indentation, allowing for 
simple renormalization of the force-indentation response to eliminate dependence 
on E. The same cannot be said for the scaling with capsid thickness t. Because 
nonlinear geometric effects are included in the finite element model, it is difficult to 
predict analytically the precise scaling with thickness. Roughly, thickness scaling 
should follow that of linearized thin-shell theory wherein force goes as the square 
of thickness F ~ t 111) , but this is only an approximation. Fig.[8]clearly shows 



the change in behavior as the thickness of the capsid is varied, while holding the 
average radius constant at 11.8 nm. 

For deformation of a thin shell, linearized elasticity theory gives a force dis- 
placement scaling F oc Et 2 d/R, where t is the shell thickness, R is the size of 
the shell, and E is the Young's modulus 111) . In general the proportionality 
factor is a function of the loading and boundary conditions of the problem. For the 
present application, Fig. |8lb) plots the the contact force, normalized by Et 2 , ver- 
sus indentation, normalized by CCMV's average radius R. Even though this model 
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Figure 8: Contact force response of CCMV over a range of capsid thicknesses. 
(Top) Force vs. AFM plate displacement for E = 250 MPa and R aV g = 11.8 nm. 
(Bottom) Dimensionless renormalization of force response according to expected 
scaling from linearized thin-shell elasticity. 
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considers geometric nonlinearities, the curves for different thicknesses nearly col- 
lapse to the linearized scaling, with an average proportionality factor of 0.8, i.e., 

F « 0.8^d. 
R 

Assuming a nominal thickness of t = 2.8nm and average radius of R = 11.8nm 
for CCMV, the Young's modulus can be estimated by relating the experimentally 
measured spring constant of the capsid to the approximate theoretical spring con- 
stant of 0.8Et 2 /R. In this way, the Young's modulus is calculated to be 280 MPa 
for the empty wild-type capsids and 360 MPa for the empty subE mutant capsids. 
Notably, these modulus values are similar to those of soft plastics, such as Teflon, 
and are slightly smaller than characteristic values for single proteins as measured 
in single molecule experiments (20). This is consistent with the notion that much 
of the local deformation of actual capsids may be sustained in the regions between 
protein subunits, such that the overall stiffness of the capsid structure is less than 
that of the individual proteins. 

The renormalized plot in Fig. |8jb) also reveals more clearly how deviations 
from linearity in the force response depend on shell thickness. For larger thick- 
nesses, the stiffening Hertz-like response near the origin is more pronounced, while 
for smaller thicknesses softening of the force response at larger indentations is 
more prominent. The mechanisms for both of these nonlinear effects are related 
to changes in the geometry of contact between the AFM tip and the capsid shell, 
which are primarily controlled by the ease with which the shell can be bent to 
accommodate the rounded tip. Notably, for thicknesses near the average physical 
thickness of the CCMV capsid (2.8 nm), these two nonlinear effects (one stiffening 
and the other softening) seem to balance each other out resulting in a nearly linear 
response. It is perhaps reasonable to conjecture that these mechanisms have roles 
in producing the linearity observed experimentally for the actual capsid. 

Nanoindentation of 029 

A model assembly similar to that applied to CCMV above has been employed to 
simulate indentation of 029, mimicking the experiments of Ivanovska et al. (|6|). As 
the same experimental setup was involved in experiments on both capsids, here the 
AFM tip geometry is again modeled as spherical with a radius of 14 nm. However, 
029 has a significantly larger capsid than CCMV, by a factor of two in diameter. 
Hence, the same 14-nm AFM tip is smaller relative to 029. In addition, the nominal 
thickness of the 029 capsid relative to its average radius is much smaller than that 
of CCMV (by a factor of four). Recalling that softening of the force response 
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Figure 9: Contact force response of 029 over a range of thicknesses, showing 
the dimensionless renormalization of the force according to expected scaling from 
linear elasticity. The average of the two principal radii (R = 23.2 nm) is used to 
scale the plate displacement. 



for CCMV was accentuated both by smaller tip size and smaller shell thickness, 
it is not surprising to see in Fig. [5] that the simulated force-indentation response 
of 029 exhibits more significant softening at moderate to large indentations over 
a range of modeled thicknesses. The response of 029 adheres reasonably well to 
the linearized scaling, however the shell-thickness-related softening nonlinearity 
appears more pronounced than for CCMV. 

From Fig. ^ the nominal capsid thickness £ = 1.6 nm, and a nominal radius 
(chosen to be the average value between the smaller and larger radii, R = 23.2 nm) 
are chosen such that an estimate of the linearized proportionality factor, and capsid 
Young's modulus can be made. From Fig.|9l for a wide range of thicknesses, a good 
estimate of the average proportionality factor is 0.6. Thus, from the experimentally 
determined spring constant of the 029 capsid, the Young's modulus is estimated to 
be E = 4.5 GPa. 

As seen in the sequence of deformed shapes in Fig. ^0] the progression of the 
von Mises stress is localized to the upper and lower regions of the capsid which 
undergo the highest deformation. However, the upper region of the capsid devel- 
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Figure 10: Indentation of ellipsocylindrical model of 029 with dimensions from 
Fig. [Q E = 1.5 GPa, and v = 0.4. Deformed capsid shapes are shown for several 
values of substrate displacement. Color contours indicate the von Mises stress. 

ops stress much more quickly, and with a higher magnitude. Due to the lack of 
axisymmetry of the capsid about the indentation axis, there is a region away from 
the region of contact with the AFM tip where the stress begins to be prominent at 
d ~ 17 nm: the inner surface of the upper portion of the capsid nearest to the cen- 
ter. However, this only occurs at displacements in the simulation that exceed the 
displacements at which failure occured experimentally. At the point when drops in 
the force are seen experimentally, d « 12 nm, the highest stresses in the simulated 
capsid are sa 760 MPa for E = 4.5 GPa. This is on par with the ultimate strength 
of typical protiens (20), suggesting that breakage of the capsid could very well be 
the cause of failure in experiments. The highest strain in the capsid is sa 11% at 
d m 12 nm, and occurs on the inner surface of the capsid directly under the AFM 
tip. However, unlike the CCMV simulations, no separation of the capsid away 
from the AFM tip is observed. 
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Discussion 



A core objective of this work has been to illuminate some of the mechanisms re- 
sponsible for the linear force-indentation responses measured for the capsids of 
viruses CCMV and 029. The acknowledgment that these nanoscale structures 
act as elastic shells prompts the question of how similar they are mechanically 
to macroscopic shells, which are modeled well by continuum elasticity theory. 
It is not clear how to answer this question directly given that a consise quanti- 
tative theoretical description of the local mechanics of the individual constituent 
capsid proteins and their interactions is currently out of reach. In the absence of 
such detailed constitutive understanding, experminents such as those performed 
with atomic force microscopy (|y,|7J) provide a coarse-grained probe to measure the 
global structural response of capsid shells. The continuum modeling in this work 
serves as a convenient coarse-grained theoretical complement, partially clarifying 
what features of the observed experimental results are generic for shells, and what 
features require detailed constitutive knowledge for explanation. 

The results here demonstrate that the generally nonlinear shape of the force- 
indentation response for a homogeneous spherical shell is affected significantly by 
the dimensions of the indentor and the shell itself. In particular, stiffening at low 
force/indentation, similar to Hertz response for solid bodies in contact, is character- 
istic of the indentation response for thicker shells with dimensions like the CCMV 
capsid. Small AFM tips and small capsid thicknesses both facilitate bending de- 
formation, which generally softens structure at higher indentations. Large tip sizes 
and shell thicknesses can have the opposite effect, leading to further stiffening of 
the force response at larger indentations. 

Consideration of nonlinearites in modeling capsid nanoindentation has an im- 
pact on the interpretation of experiements. Previous linearized elasticity models 
constructed by Ivanovska et al. (Jg) estimated a Young's modulus for 029 of 1.8 
GPa. This result appears mainly to be a function of two things: a proportionality 
factor found analytically for a spherical shell with point loads applied at the top and 
bottom, and a finite element model also designed with point loads with the capsid 
modeled as a geodesic ellipsoid. Together, these factors lead to an estimate of the 
Young's modulus which is smaller than that of the present model. Clearly such 
difference is to be expected with the point loading exchanged for a finite-sized tip. 
However, both the previously reported value and the value reported here show the 
029 capsid to have a stiffness in the range of 1 — 5 GPa. 

The continuum shell models presented here may also help to explain how the 
linearity of the experimentally measured force-indentation responses of 029 (@) 
and CCMV (7) can persist despite the pervasiveness of nonlinearity in the system. 
The dimensions of the AFM tips and capsids in recent experiements on CCMV are 
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in a range in which the nonlinear stiffening and softening effects seem to balance 
each other so that the modeled elastic force-indentation response is close to lin- 
ear. This result suggests that similar experiments done with AFM tips or capsids 
of other dimensions might reveal force-indentation responses which are nonlinear 
even in the elastic regime. However, though the measured response of 029, which 
has significantly different dimensions from CCMV, is also linear, the continuum 
models here still show noticeable nonlinearity. This indicates a deficiency of the 
present continuum models for precise prediction of capsid behavior. This could 
perhaps be remedied by a more accurate representation of the geometry of the 029. 
Alternatively, the models could be improved with a more informed choice of con- 
stitutive law, although the constitutive models employed here, which are capable 
of significant nonlinearities, did not significantly affect capsid shell response. 

Another drawback to the isotropic, homogeneous, geometrically simple model 
is that it does not produce drops in contact force that are experimentally observed 
at sa 20 — 30% deformation. These drops are signs of failure of capsid structure, 
which in general could be triggered by geometric instability (buckling) or local 
failure of the material (e.g., bond-breaking). Over the relevant range of dimen- 
sions, the models in this work did not exhibit any geometric instabilities. How- 
ever, the present model does not take into account any of the geometric complex- 
ities present in the molecular structure; although spherical in an average sense, 
the CCMV capsid is faceted into distinct regions of pentamers and hexamers and 
there are specific and finite points where the subunits are joined through chemi- 
cal bonding. The same is true for 029. These geometric inhomogeneities could 
change the buckling characteristics of the structure. Also, the present model ne- 
glects the effects of "pre-stressing" consistent with the hypothesis of Lidmar et al. 
(|2|) that the icosahedral vertices of spherical capsids act as disclinations in an other- 
wise hexagonal protein lattice. In fact, based on this hypothesis, the recent coarse- 
grained molecular-dynamics simulations of Vliegenthart and Gompper (23) reveal 
that buckling, marked by force drops, is indeed possible during capsid indentation. 

Because the nonlinear elasticity model did not turn out to predict any drop 
in force due to geometric instability, it may be inferred that the mechanism of 
failure during the experiment is likely to be, at least in part, capsid breakage. To 
test this hypothesis, a more complex model would be needed that describes the 
physics of interactions between the subunits and the more variated structure of the 
capsid. A plausible scenario is that the capsid structure fails at the noncovalent 
bonds between subunits. It is possible that a true atomic model would be required 
to capture points of failure nucleation/initiation. 
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